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Heterogeneity of cells may explain allometric scaling of metabolic rate 
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Abstract 

The origin of allometric scaling of metabolic rate is a long-standing question in biology. Several models have 
been proposed for explaining the origin; however, they have advantages and disadvantages. In particular, 
previous models only demonstrate either two important observations for the allometric scaling: the variability 
of scaling exponents and predominance of 3/4-power law. Thus, these models have a dispute over their 
validity. In this study, we propose a simple geometry model, and show that a hypothesis that total surface 
area of cells determines metabolic rate can reproduce these two observations by combining two concepts: 
the impact of cell sizes on metabolic rate and fractal-like (hierarchical) organization. The proposed model 
both theoretically and numerically demonstrates the approximately 3/4-power law although several different 
biological strategies are considered. The model validity is confirmed using empirical data. Furthermore, the 
model suggests the importance of heterogeneity of cell size for the emergence of the allometric scaling. The 
proposed model provides intuitive and unique insights into the origin of allometric scaling laws in biology, 
despite several limitations of the model. 

Keywords: Allometric scaling. Metabolic rate. Geometry model. Fractal-like organization. Log-normal 
distribution 


1. Introduction 


M etabolic processes are essentia l for physiological functions and responsible for maintaining life (|Takemoto 


20121 : iTakemoto and OosawaLl2012f) . The relationship between metabolic rate B and body mass M is an im¬ 
portant and interesting topic of scientific inquiry not only for researchers in the field of basic biology but also 
for investigators in ecology and medic al research, and it is well known to approximately obey a power law 
( West et ah . 2002 : Brown et ah . 200^: B oc . This allometry is positioned as a significant equation in 
both biol ogy and ec o logy b ecause it is useful for un derstanding and f or estimating the energy metabolism, 
lifespan ( Speakman . l2005ll . and animal space use ( Jetz et ah . l2004h : in particular, determination of the 
scaling exponent 7 is a long-standing question. 


A pioneering study includes the surface law found by Rubner in the 1880s ((White and Seymoud . 1200311 : 
in particular, Rubner reported that metabolic rate is proportional to in mammals. The Rubner’s 

surface law is immediately deriyed when assuming simple geometric and physical principles: metabolic rate 
(e.g., heat production rate) is proportional to the rate of energy (e.g., heat) dissipated through body surface 
because of homeostasis. 

Contrary to the prediction from the surface law, in 1932, Kleiber proposed t hat metabolic rate is proper- 


Uontrary to tne prediction trom tne surtace law, m t vleiber proposed mat metabolic rate is propor¬ 
tional to ( White and Seymour . 20031:IWest et al.l . 120021 : [^eakmanl . 2005 : Brown et al.l . 200411 . Further 


studies have confir med that the Kleib er’s (i.e., 3/4-power) law is predominant at least in plants and ani¬ 
mals (reviewed in ( Savage et al.l . 2004)): this allometric scaling is observed in wide-ranging organisms (i.e., 
microorganisms to elephants). 
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West et alJ (|l997fl have proposed a model (West-Brown-Enquist (WBE) model) for explaining the origin 
of the 3/4-power law. This model assumes that oxygen and nutrients are transported through space-hlling 
fractal networks of branching tubes, in which the number of capillaries (leaves in the case of plants) is 
proportional to metabolic rate. Since the WBE model clearly illustrates the Kleiber’s law, it is frequently 
used for understanding the allometric scaling. 

In a ddition to th i s, fur ther studies have proposed sever al alternative m o dels based on the f ourth dimension 
of life ( West et"^. 1999ll. optimal transpor t networks ( Banavar et ah . 19991 20101 2014 1. and quantum 
metabolism ( Demetrius and Tuszvnskil 201011 . 

Price_eLaD (2012 ) have ca refully evaluated the 


However, these models have several limitatio ns (e.g 


WBE model using empirical data). In particular, Kozlowski and Konarzews5 ( 2004h have questioned the 
universality of the 3/4-power law. In f act, s everal studies suggest that the scaling exponent 7 is vari¬ 
able. For example. IWhite and Sevmouii (2003) have reported that metabolic rate is proportional to 
in mammals when considering body temperature, digestive state, and phylogeny. IParveau et al.l (|2002ll 
have illustrated that the approximately 3/4-power law results from the sum of the influences of multiple 
contributors to metab olism and control (i.e., the sum of allometric relationships observed in a number of 
biological processes). iReich et al.l ( 20061) found that a linear relationsh ip between the rate of respiratory 
metabolism and body mass (i.e., B oc M) is observed in plants although lEnauist et al.l (120071 1 have refuted 
this co nclusion. Similarly, the observed similar mean mass-specific metabolic rates across life’s major do¬ 


mains (|Makarieva et al.l . 1200811 also implies that B (x M. To explain the variability of th e scaling exponent 


however, the WBE model can be modified ( Price et al. . 20071 : Kolokotrones et HI . I2OIOII . 


_ O f particular intere s t is th e fact th at the scaling exponent varies according to cell size (|Kozlowski et al.l . 

2003 : Starostova et al.L l2009l) . Thus, iKozlowski et al.l ( 2003ll have proposed a simple model (Kozlowski- 
Konarzewski-Gawelczyk (KKG) model) based on cell size, inspired by the fact that a large part of standard 
meta bolic costs are spent preserving ionic gradients on cell membranes ( SzarskiL 19831 : Porter and Bran^ 
199311 . This model considers an extension of the Rubner’s surface law: a hypothesis that metabolic rates 
are determined by total surface area of cells rather than body surface. The KKG model can explain the 
variability of the scaling exponent; especially, the exponent can fall within the range between 2/3 and 1. 
For example, B oc when cell size is proportional to body size. On the other hand, B oc M when cell 

size is indepe ndent from body size. The va riability of the scaling exponent has already well known in terms 
of geometry (lOkid. I 2 OI 3 I : iHirst et'H! . 2014ll . 

iBrown et al.1 (|2005l) hav e argued against the K KG model. Since cell size is almost independent from 


body size, as explained bv IKozlowski et al.l (|2003D . the KKG model s uggests B oc M, indicating that it 


contradicts the predominance of the 3/4-power law (jSavage et al.l . 120041) . 


Is the KKG model or a hypothesis that the total surface area determines metabolic rate not really 
useful for understanding the allometric scaling law? In this study, we propose a simple geometry model, an 
extended KKG model, and show approximately 3/4-power law can be also emerged from this hypothesis by 
considering the concept of fractal-like (hierarchical) organization of the WBE model. This result suggests 
that the hypothesis is evidently useful for understanding allometric scaling of metabolic rate. Moreover, our 
model explains both validity of the scaling exponent and ubiquity of the 3/4-power laws, and it suggests 
that the allometric scaling of metabolic rate is affects by cell size distributions, rather than cell sizes. 


2. A geometry model 

The KKG model assumes that organisms consist of uniform distributed isometric cells; however, such an 
assump tion may be unsat isfied because organisms show fractal-like (hierarchical) organization, as pointed 
out bv IWest et al. (1997). Furthermor e, heterogeneous dis tributions such as log-normal distribution are 


widely observed in real-world systems (jLimpert et al.l . 1200111. Multiplicative effec ts and hierarchical orga¬ 


nizations are know n to generate log-normal distrib utions ( Kobavashi et al. . 201lh . In particular, a simple 


model of fracture ( Yamamoto and Yamazaki . 20131) . originally proposed by A. N. Kolmogorov in 1941, is 
an instructive example, and it considers hierarchical divisions of one rod. Thus, we expected that such 
heterogeneous distributions of cell sizes resulting from the hierarchical organization lead to H oc with 
7 < 1 even if (mean) size of cell sizes is almost independent from body size. 
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Inspired by this fracture model, we propose a geometry model, called fractal-like cube division (FCD) 
model (Fig. [T]). The FCD model can be interpreted as a 3-dimensional version of the fracture model; however, 
note that division processes are slightly different between the fracture model and FCD model because of the 
consideration of fractal-like (hierarchical) organization of cell sizes (e.g., tree branches are more frequently 
divided than the tree trunk is). 

We first explain a simple case of the FCD model. In the simple FCD model, the cube with the length of 
L (Fig. [T]4) is divided according the following procedures, (i) The cube is divided into 8 cubes (Figs. [TJ3 
and 121^). (ii) A randomly selected cube is further separated into 8 cubes (Fig. [IP), (iii) The procedure (ii) 
is repeated until t = T (Fig. HP). 



? = 0 r=l r = 2 r = 7’ 


Figure 1: Schematic diagram depicting the fractal-like cube division (FCD) model. (A) Cube with the length of L. (B) The 
cube of (A) is divided into 8 cubes (cuboids) with the length of L/2. A cube (cuboid) selected at random (i.e., with the 
probability of 1/8) is indicated by gray. (C) As in (B), the selected cube (cuboid) is separated into 8 cubes (cuboids). (D) A 
state of the cube of (A) after T divisions. 


However, this division rule of the simple FCD model may be by-the-numbers. To consider more flexible 
divisions, we next propose a general FCD model. This model considers a division of the cuboid Ch,w,d 
composed of the height of /i, width of rc, and depth of d, in the procedure (ii). In particular, this cuboid 
is divided into 8 cuboids according to the parameter p drawn from a probability distribution with the 
range from 0 to 1 (e.g., uniform distribution C/(0,1)) (Fig. HP): Cph,pw,pd, Cph,pw,ii-p)d, Cph,(i-p)w,pd, 

^ph^{l—p)w,{l—p)dj C{l—p)h,pw,pdi C{l—p)h,pw,{l—p)d^ ^{l—p)h,(l—p)w,pd^ and p)c/- 

The general FCD model is equivalent to the simple FCD model when p = 0.5 at all times. 



Figure 2: Division rules of the simple FCD model (A) and general FCD model (B). 


3. Results and discussion 

3.1. Analytical solutions of geometric parameters 

We here provide analytical solutions of geometric properties of the simple FCD model; however, the 
analytical results may be applicable to the general FCD model (see Sec. 13.211 . 
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To obtain an analytical solution of the surface area of the simple FCD model using a mean field ap¬ 
proximation, we consider the time evolution of the number nD{t) of cubes with the length of (1/2)^L when 
D > 0 (i.e., t > 0), where D denotes the cube division number. 

no{t) increases by 8 when a cube with the length of is selected with the probability of 

nD-i{t)/N{t), and it decreases by 1 with the probability of nD{t)/N{t), where N{t) denotes the total 
number of cubes at time t. Thus, noit) is described as: 




npjt) 
Nit) ’ 


( 1 ) 


where A^(t) = 7t -I- 1 ~ 7t (when t0). 

Since no(t) = 0 when t > 0, the time evolution of ni(t) is approximately written as 

d , , n-iit) 

The solution of Eq. m with the initial condition of ni(l) = 8 is 

We solve the recursive equation Eq. © with this initial condition Eq. m, and finally obtain 


( 2 ) 

( 3 ) 


noit) = 


8 


D 


(Int) 


D-l 


7D-HI/7 

Therefore, the surface area A of the simple FCD model at time T is described as 

2 


Dr, 


D 


nniT) 


A= ^ 6 

D^l 

_ ior2rTxl/7 ^(-^max) f lu T) 

r(D^ax) ’ 


( 4 ) 


( 5 ) 


where Dmax denotes the maximum number of cube division (i.e., the minimum cell length is (1/2)^”“). 
r(x) and r(a;, j/) indicate the Euler’s gamma function and incomplete gamma function, respectively. 

Assuming D^ax ^ 0 (i.e., the minimum cell length (1/2)^““ -C L), the above equation is approximately 
written as 

A = (6) 

The relative frequency P{D) (i.e., probability distribution) of cubes with length of [1/2)^L (at time T) 
is approximately equivalent to npiT)/[IT)-, thus, we have 


PiD) ^ 

where A = (8/7) InT. Equation (O indicates that PiD) is a Poisson distribution with the mean D of A -I- 1, 
suggesting a geometric parameter x of divided cubes (e.g., length I, surface area a, volume v) follows a 
log-Poisson distribution because D oc lux. 


3.2. Numerical evaluations 


To evaluate the validity of these theoretical results, we performed numerical simulations. 

We first compared the relationship between the surface area A and volume V (i.e., = L^) between our 
theory and numerical simulation. Supposing that T = , we immediately obtain the relationship between 

A and V from Eq. 


A = 12V 


( 8 ) 
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The prediction from Eq. (I5|) is in excellent agreement with the numerical results in both cases of the 
simple and general FCD models (Fig. [3]). In this case, the parameter p in the general FCD model was drawn 
from the uniform distribution C/(0,1) when dividing a cuboid. 

We also confirmed that Eq. ([8|) can estimate the A-V relationships numerically obtained from the general 
FCD model, in which the parameter p is sampled from the beta distribution Beta(2,2), a bell-shaped 
distribution whose random variables range between 0 and 1, when separating a cuboid. The prediction 
accuracy in this case is almost equal to that in the case of the uniform distribution. 

These results indicate that the analytical solutions of the simple FCD model can be applicable to the 
general FCD model when the parameter p is a random variable drawn from a probability distribution with 
the mean of 0.5 and limited variance. The applicability of the analytical result obtained from the simple 
model to the general model is due to equivalents of the central limit theory. 



Volume 


Figure 3: (Color online) Relationships between surface area and volume of the simple FCD model (A) and general FCD model 
(B), in which the parameter p is drawn from the uniform distribution 17(0,1). These simulations consider that T = [V^J. 
Symbols denote numerical results, which are averaged over 100 realizations. Lines are drawn from Eq. lEJi. Insets represent 
the same plots in normal scale. 
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We next compared the predicted relative frequency (probability density) of geometric parameters with 
numerical simulations. As an example of geometric parameters, we focused on cube volumes v. Using the 
relational expression D = — ln(u/U)/ln8, the relative frequency can be estimated from Eq. ([7]). 

As shown in Fig. |4j our theoretical predictions are in good agreement with numerical results although 
the prediction accuracy becomes lower for abundant cubes with a smaller volume due to the several approx¬ 
imations (see Sec. 13.11) . However, this limitation poses little problems when estimating the total surface 
area A of the model because the theoretical results are in excellent agreement with numerical results (Fig. 

El). 

We also confirmed that Eq. © can predict numerical results in both cases of cube lengths I and cube 
surface areas a. Substituting the relational expression D = — ln(l/L)/ln2 {D = — ln(a/Ao)/ln4, where 
Aq = 6 L^) into Eq. (I?]), we can obtain the theoretical distributions of I (a). 



Figure 4: (Color online) Relative frequencies of cube volumes v in the simple FCD model with L = 40 and T = 100, 1000, and 
10000. Symbols denote numerical results, which are averaged over 100 realizations. Error bars indicate the standard deviations. 
Lines are drawn from Eq. 0 using D = — ln(u/y)/ InS. 


3.3. Comparison with empirical data 

Assuming B oc A and M oc V ( Kozlowski et ah . 20031: White and Seymour . 2003), we can discuss the 
origin of the allometric scaling of metabolic rate using the FCD model; however, the parameter T needs to 
be estimated by considering some constraints. We here estimate T using empirical results, and derive the 
allometric scaling laws of metabolic rate. 


For example, two extremely possible strategies ([Savage et al 


20071) may be useful: (i) average cell (cube 
or cuboid) mass (m) remains constant while average cellular metabolic rate ( 6 ) changes, or (i i) {b) remains 


fixed while (m) varies, where {x) indicates the mean of x. According to a previous study ( Savage et al.l . 
20071) . we assume that (m) oc (u) and ( 6 ) oc (a), where {v) and (a) correspond to average cellular volume 
and average cellular surface area (a), respectively. In this case, the parameter T can be estimated as follows: 
In particular, the strategy (i) leads to T cx M because (m) oc V/T. Moreover, the strategy (ii) derives 
T oc because (a) = l{7T). Therefore, we find 


B oc 


^17/21 ^ jy-0.81 
^ 7/9 ^ ^ 0.78 


(strategy (i)) 
(strategy (ii)) 


(9) 


On the other hand, if T is proportional to the number of cells Nc because of the definition, our model 
predicts Nc oc M and W oc in the case of the strategies (i) and (ii), respectiv ely. Theoretical 

predictions are in agreement with the observed results (see Table 1 in ( Savage et al.l . l2007l) ). 

However, since the distributions of a and v are log-normal like (Fig. (Hand Eq. ©), it remains possible 
that the estimation of T using means is unsuitable. When considering the multiplicative process such as 
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the FCD model (cell divisions), from a statistical perspective, geometric means should also be considered. 
Thus, we further considered 2 strategies by modifying the strategies (i) and (ii): (iii) logarithmic mean of 
cell mass (Inm) remains constant while logarithmic mean of cellular metabolic rate (Inti) changes, or (iv) 
(Inti) remains fixed while (Inm) varies. 

(In v) and (In a) in the simple FCD model are obtained as 


and 


(Inn) = In 
£ 1=1 



= 31 In - - -lnTln2 I 

V 2 7 J 


(Ino) = In 
£ 1=1 
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= 2 


/ T O 

In-lnTln2 

V 2 7 



( 10 ) 


( 11 ) 


respectively. 

Supposing that (In b) oc (In a) and (In m) cx (In v ), both strategies (iii) and (iv) lead to T oc . 

Thus, we obtain 

B oc 1173+24^2 = (strategies (iii) & (iv)). (12) 

The scaling exponent ranges from 0.73 to 0.81 according to the biological strategies, and it may fall with in 


95% confidential intervals of scaling exponents estimated from empirical dat a (e.g., see ( Savage et ah, 2004ll b 


In addition to this, this result is consistent with a conclusion, derived bv iKolokotrones et al.l ( 2010ll . that 


the scaling exponent is around 0.8 rather than 3/4 when only focusing on mammals with higher M. This 
finding may not contradict the WEB model because t his model is known to show the scaling exponent of 
0.81 (not 3/4) when considering finite-size corrections ( Savage et al.l . 12008 1. 


The relationship between the C-value (i.e., genome size) and M I Kozlowski et al.l . l2003ll is also useful 
for estimating the paramete r T ; in p a rticul ar, the C-value has been known to show an approximately linear 


relationship with cell size (Gregory, 200 lll . This previous study has reported that C-value oc with 


P G (—0.05,0.06); thus, it encourages an extension of the strategies (i) and (iii). In particular, we considered 
the following extended strategies: (I) (m) oc and (III) (Inm) cx /31nM. When /3 = 0, the strategies (I) 
and (III) are equivalent to the strategies (i) and (iii), respectively. 

Using the relational expression (m) (xVjT and Eq. (flllll . we can obtain the allometric relationships of 
B based on the strategies (I) and (III), respectively. In particular, the scaling exponents are described as 


and 


respectively. 


7 = 


^=3 


17-3/3 

~21 

1-/3 

24 In 2 


(strategy (I)) 


(strategy (III)), 


(13) 


(14) 


We evaluated these theoretical predictions with the empirical data ( Kozlowski et al. . 2003ll (Fig. [S]). 


For comparison, in addition, we also investigated the 7-/3 relationship predicted from the KKG model. 
The total surface area of the KKG model is described as cx (m)^/^/Vc, where Nc oc Mjpm). Since 
ln(m) = (Inm) cx /31nm because all cell sizes are identical (i.e., the variance of cell size is zero) in the KKG 
model, the 7-/3 relationship, derived according to the strategies (I) and (III), is the same as follows: 


7 = 1 — ^ (KKG model). 


(15) 









































Our theoretical predictions are in good agreement with the empirical data, and they show higher predic¬ 
tion accuracy, compared with the KKG model. Our theory can also explain a negative correlation between 
metabolic rate exponent 7 and C-value exponent j3. However, the theoretical prediction of this tendency 
may be qualitative in this case. In particular, the FCD model could not quantitatively predict the 7-/3 
relationships for a larger 13. This may be because the exponent (3 was estimated from average cell size (i.e., 
(m)) calculated using a few number of cells. That is, it remains possible that the assumption that (m) oc 
and (Inm) oc (3\nM is not satisfied. A more careful examination of the theoretical predictions may require 
an estimation of /3 using a larger number of cells in each organism. 
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Figure 5: (Color online) Relationship betwe en metabolic rate (MR) exponent 7 and C-value exponent /3. 
empirical data, obtained from the study bv iKoztowski et al.l (l2QQ3lh and they indicate lineages (orders), 
standard errors. The solid lines are drawn from Eqs. and Gl (our (FCD) model), respectively, 

obtained from Eq. 1I15II (KKG model). 


Symbols denote the 
Error bars indicate 
The dashed line is 


We can also consider an extension of the strategies (ii) and (iv) according to the M-dependence of cell 
size ( Kozlowski et ah . 200,3h : (II) (a) oc and (IV) (Ino) oc /31nM. Note that /3 does not correspond to 
the C-value exponent in this case. When /3 = 0, the strategies (II) and (IV) are equivalent to the strategies 
(ii) and (iv), respectively. 

According to the relational expression (a) = 12M^/^T^/’^/(7T) and Eq. (fTTl) . we find the allometric 
relationships of B based on the strategies (II) and (IV), respectively. Especially, the scaling exponents are 
obtained as 

14-3/3 , 

7 = ——— (strategy (II)) 


and 


respectively. 


7 = — 

3 


18 


2-3/3 
48 In 2 


(strategy (IV)), 


(16) 


(17) 


3 . 4 . Ubiquity of the 3/4-power law 

The strategies (I)-(IV) may be mixed among real-world living organisms. Furthermore, a huge variety of 
(3 may be observed. What kind of distributions does the metabolic rate exponent 7 follow in this situation? 

To consider various possible situations, we assumed the following condition: the strategies (I)-(IV) are 
uniformly observed, and /3 is drawn from the uniform distribution [/(—1,1). 









































The 40000 trials performed under this condition showed a bell-shape distribution of 7 with the mean of 
0.76 (the median of 0.75) (Fig. [5] A) although the strategies and /3 are uniformly selected. 

In terms of evolutionary dynamics, it may be more realistic to consider that /3 follows a Gaussian-like 
distribution. We also performed a similar analysis by assuming that /3 is sampled from a normal distribution 
with the mean of 0 and variance of 1 (i.e., A^( 0 , 1 )), and again found that the metabolic rate exponent follows 
a bell-shape distribution of 7 with the mean of 0.76 (the median of 0.75) (Fig. IHlB). 

This result implies that the 3/4-power law can be emerged from a number o f different biological strate¬ 
gies. Although a previous geometric perspective ( OkieL [201^ iHirst et al.1 . l2014h shows that the metabolic 
rate exponent 7 can be variable, it cannot answer why the 3/4-power law is universally found in nature 
( Savage et al.l . l200'3l . The FCD model may answers this question. 



0.60 0.65 0.70 0.75 0.80 0.85 0.90 0.95 
MR exponent 



MR exponent 


Figure 6 : (Color online) Distribution of metabolic rate (MR) exponent 7 obtained from the FCD model according the strategies 
(I)-(IV), in which f) is drawn from the uniform distribution U{—1, 1) (A) and the standard normal distribution N{0, 1) (B). 


4. Conclusion 

We proposed a simple geometry model for the emergence of allometric scaling of metabolic rate. Several 
analytical solutions of this model of geometric parameters were provided (see Sec. ED, and their validity 
was confirmed by numerical simulations (see Sec. 13.21) . 

Using this model, we discussed a possible origin of allometric scaling of metabolic rate, and concluded 
that the approximately 3/4-power law can be also derived from the context of total surface area of cell sizes 
by considering fractal-like (or hierarchical) organization of living organisms, contrary to criticisms (e.g.. 
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(|Brown et alJ . l2005ri l. The theoretical predictions obtained from this model are in good agreement with 


empirical results (see Sec. Unlike previous models, in particular, this model can explain the important 
observation s: the metabolic rate exponen t 7 shows a central tendency of 3/4 (see Sec. 13.41) although it can 
be variable ( Okie! . 2013 : Hirst et ^ . 12014 1. These results implies that a hypothesis total surface area of cells 
determines metabolic rate evidently shows promise for explaining the origin of the allom etric scaling. 

The FCD mod el was proposed inspired by the KKG model ( Kozlowski et al.l . 12003 ) and WBE model 


(West et al 


19971) : in particular, it considers the effect of cell size (KKG model) on metabolic rate and the 
concep t of f ractal-like organizat ion (WBE model) (see Sec. [2]for details). Although lKozlowski and Konarzewski 
( 2004 ) and Brown et al.l ( 20051) debated over the validity of each other’s models, the EDC model suggests 
that these models are complementary rather than conflicting. In particular, the EGD model shows a fractal- 
like organization (heterogeneous distribution) of cells is required that when deriving the approximately 
3/4-power law in the context of geometry (the KKG model). 

The EGD model implies the importance of heterogeneous distributions (e.g., log-normal distributions of 
cell sizes for the allometric scaling. In this paper, we could not test this hypothesis because c ell size distribu¬ 


tions are less understood. However, log-normal distributions are widely observed in nature (jLimpert et al 
I 2 OOII) . and several previous studies suggest that the size o f red cell, which is rel ated to metabolic rate, 
approximately follows a log-normal distribution ( Bell . 1985t McLaren et al.l . [T99lh . Therefore, we expect 
that this hypothesis provides a new perspective for the origin of metabolic scaling. Although several studies 
( Kozlowski et al.l . 20031 : IStarostova et al.l . l2009l) have reported the importance of cell size for scaling laws of 
metabolic rate, there is little study focusing on such size distributions. Our study encourages an evaluation 
of the relationship between the scaling laws of metabolic rate and cell size distributions (i.e., single-cell 
measurements). 

However, more careful examinations may be required. The relationship between the surface area and 
metabolic rate is still unclear because of several reasons. Eor example, a number of factors influence metabolic 
rates. The rate of supply and transportation of resources within the body is considered as the fundamental 
constraint on metabolic rate; however, the EGD model does not consider this constraint. Although the total 
surface area of mitochondrial membranes on which the majority of ATP synthesis occurs may be correlated 
with metabolic rate, it is not conformed whether such a correlation is also observed in the other types of 
cells. 

Even though limitations are inherent in our model due to the simplicity of the model, we believe that our 
model still serves to provide intuitive and unique insights into the origin of allometric scaling in metabolic 
rate. 
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